Geometry of the Computational Singular Perturbation Method
نویسندگان
چکیده
The Computational Singular Perturbation (CSP) method, developed by Lam and Goussis [Twenty-Second Symposium (International) on Combustion, The Combustion Institute, Pittsburgh, 1988, pp. 931–941], is a commonly-used method for finding approximations of slow manifolds in systems of ordinary differential equations (ODEs) with multiple time scales. The validity of the CSP method was established for fast–slow systems with a small parameter ε by the authors in [Journal of Nonlinear Science, 14 (2004), 59–91]. In this article, we consider a more general class of ODEs which lack an explicit small parameter ε, but where fast and slow variables are nevertheless separated by a spectral gap. First, we show that certain key quantities used in the CSP method are tensorial and thus invariant under coordinate changes in the state space. Second, we characterize the slow manifold in terms of these key quantities and explain how these characterizations are related to the invariance equation. The implementation of the CSP method can be either as a one-step or as a two-step procedure. The one-step CSP method aims to approximate the slow manifold; the two-step CSP method goes one step further and aims to decouple the fast and slow variables at each point in the state space. We show that, in either case, the operations of changing coordinates and performing one iteration of the CSP method commute. We use the commutativity property to give a new, concise proof of the validity of the CSP method for fast–slow systems and illustrate with an example due to Davis and Skodje. Dedicated to Professor Alexander Gorban on the occasion of his sixtieth birthday
منابع مشابه
Near-Optimal Controls of a Fuel Cell Coupled with Reformer using Singular Perturbation methods
A singularly perturbed model is proposed for a system comprised of a PEM Fuel Cell(PEM-FC) with Natural Gas Hydrogen Reformer (NG-HR). This eighteenth order system is decomposedinto slow and fast lower order subsystems using singular perturbation techniques that provides tools forseparation and order reduction. Then, three different types of controllers, namely an optimal full-order,a near-opti...
متن کاملA Discrete Singular Convolution Method for the Seepage Analysis in Porous Media with Irregular Geometry
A novel discrete singular convolution (DSC) formulation is presented for the seepage analysis in irregular geometric porous media. The DSC is a new promising numerical approach which has been recently applied to solve several engineering problems. For a medium with regular geometry, realizing of the DSC for the seepage analysis is straight forward. But DSC implementation for a medium with ir...
متن کاملDetecting the location of the boundary layers in singular perturbation problems with general linear non-local boundary conditions
Singular perturbation problems have been studied by many mathematicians. Since the approximate solutions of these problems are as the sum of internal solution (boundary layer area) and external ones, the formation or non-formation of boundary layers should be specified. This paper, investigates this issue for a singular perturbation problem including a first order differential equation with gen...
متن کاملOn Approximate Stationary Radial Solutions for a Class of Boundary Value Problems Arising in Epitaxial Growth Theory
In this paper, we consider a non-self-adjoint, singular, nonlinear fourth order boundary value problem which arises in the theory of epitaxial growth. It is possible to reduce the fourth order equation to a singular boundary value problem of second order given by w''-1/r w'=w^2/(2r^2 )+1/2 λ r^2. The problem depends on the parameter λ and admits multiple solutions. Therefore, it is difficult to...
متن کاملInstitute for Numerical and Computational Analysis Dublin, Ireland An adaptive uniformly convergent numerical method for a semilinear singular perturbation problem
A singularly perturbed semilinear two point boundary value problem is considered, without any restriction on its turning points. A difference scheme is presented for solving this problem on an arbitrary locally quasiuniform mesh. It is shown that the solution of the scheme is first order accurate, uniformly in the perturbation parameter, in a discrete L1 norm. Numerical results are presented fo...
متن کامل